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r-> . 

It is extremely important to devise a reliable method to extract spectroscopic factors from transfer 
cross sections. We analyse the standard DWBA procedure and combine it with the asymptotic 
normalisation coefficient, extracted from an independent data set. We find that the single particle 
parameters used in the past generate inconsistent asymptotic normalization coefficients. In order 
to obtain a consistent spectroscopic factor, non-standard parameters for the single particle overlap 
functions can be used but, as a consequence, often reduced spectroscopic strengths emerge. Different 
■»— 4- ■ choices of optical potentials and higher order effects in the reaction model are also studied. Our 

(Sj ; test cases consist of: 14 C(d,p) 15 C(g.s.) at E l d ab = 14 MeV, 16 0(d,p) 17 0(g.s.) at E l d ab = 15 MeV and 

40 Ca(d,p) 41 Ca(g.s.) at E l d ab — 11 MeV. We underline the importance of performing experiments 
H ' specifically designed to extract ANCs for these systems. 
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t^. . I. INTRODUCTION 

O ; 

The shell model formalism first introduced Spectroscopic Factors (SF) to describe the shell occupancy In 
particular, Sfj 1 , the single particle SF, is defined as the norm of the overlap function of a nucleus (A+l) in a 
particular state i with a nucleus A in a state /, where the valence nucleon is in an orbital with orbital and total 
angular momentum These SFs have been extensively compared with those extracted from reactions. At present 

ab-initio calculations are improving the accuracy of the calculated SFs (e.g. [1, One would eventually like to have 
a very accurate probe that could test the predictions of these models and could disentangle the relevant elements of 
•rH . the NN force that are still missing, especially when moving toward the driplines. 

From the experimentalist point of view, a spectroscopic factor is a ratio of measured to predicted differential 
cross sections. Phenomenological spectroscopic factors are extensively used in a variety of topics, from nuclear 
reactions to astrophysics or applied physics, yet the procedure for their extraction from the data has remained 
essentially the same for decades. Most of the work in the sixties and seventies used direct transfer reactions, such 
as (d,p), (d,t), ( 3 He, d), ( 3 He, a), as the central tool [1, H, @|. However transfer analyses have a reputation of large 
uncertainties. 

Other methods to extract spectroscopic information include the (e,e'p) reactions or the nuclear knockout in inverse 
kinematics. For stable nuclei, there have been many electron knockout experiments. They have provided SFs with 
rather small error bars [3]. The analysis of these measurements along a wide mass range, for single particle states 
with expected SFs close to unity, show an overall reduction of the SF (~ 0.6). The source of this reduction is not yet 
well understood although one expects NN short or long range correlations which are not included in the present 
day non-ab-initio shell model to contribute. 

Nuclear Knockout in inverse kinematics using radioactive beams is a new probe introduced at the NSCL [§]. A 
systematic study of spectroscopy on a variety of nuclei, ranging from the stability valley to the proton dripline, have 
been performed using this technique. It is found that the measured spectroscopic factors suffer from a reduction 
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relative to the Shell Modelpredictions and the reduction factor changes with binding energy [ltj. Again, it is not 
clear where this quenching [l| comes from. Of course, as (e,e'p), this technique is only suitable for studying the single 
particle structure of the ground states of nuclei. 

The (e,e'p) cross section is sensitive to the structure all the way to the inside of the nucleus, whereas the transfer 
is typically peripheral and surface peaked. The knockout reactions with radioactive beams have been performed in a 
kinematical regime where the eikonal approximation can be used, to simplify the reaction theory and to reduce the 
reaction model uncertainties [f}]. Nuclear knockout results |lj are in agreement with the (e,e'p) SFs for the tested 
stable closed shell nuclei [ll|. Transfer and (e,e'p) have also been compared Q for 48 Ca. Earlier transfer studies 
provided a spectroscopic factor for the last proton close to unity, but a reanalysis by including finite-range and 
non-local effects in the transfer reaction model, show that the transfer result can be brought down to 0.6, the (e,e'p) 
value. 

The spectroscopic factor is the norm of the overlap function, which peaks well inside the nuclear radius Rn- As 
the (e,e'p) reaction can probe the nuclear interior, it is suitable for extracting SFs as long as the reaction mechanism 
is well understood. However there are some problems with the high momentum transfer component, associated with 
probing the inside of the nucleus, because then the Born approximation is not valid [T3|. For exotic nuclei near or 
on the driplines, transfer reactions are a unique tool and, hence, can have a large impact in the programs of the new 
generation rare isotope laboratories. 

The standard framework for analysing transfer data with the intent of extracting SFs is the distorted-wave Born 
approximation (DWBA). Overall, it has been very successful in describing angular distributions at forward angles 
and less so for the larger angles where higher order become more important. The SF is the normalization needed for 
the calculated DWBA differential cross section to match the experimental one at forward angles (e.g. [l3l. ILH. fl5|). 
The uncertainty of the extracted SF resulting from the normalization of the DWBA cross section is assumed to be 
~ 30%, even if the statistical errors are low. The reasons for this inaccuracy are typically attributed to ambiguities 
in the optical potentials, the inadequacy of the DWBA reaction theory, or the dependence on the single-particle 
potential parameters. Recently, systematic studies on 12 C(d,p) 13 C have shown that it is possible to bring SFs into 
conformity using a global optical potential prescription [H,E3|> whereas arbitrary choices of the optical potential will 
hold disparate results. As to the reaction mechanism, there are many studies on the validity of DWBA (e.g. [18]) 
and typically the reaction mechanism needs to be checked case by case. In addition, Hartree-Fock densities have been 
suggested as a mean to constrain the single particle parameters Therein, the single particle radii for the Ca 

isotopes were adjusted to reproduce the known rms matter radii, and zero range DWBA calculations were performed 
to extract the SFs. 

In (2(J, a combined method of extracting SFs from transfer reactions was introduced. This method can also be 
applied to breakup and (e,e'p) reactions. The combined method, which is based on the introduction of the asymptotic 
normalization coefficients (ANC) into the transfer analysis, allows one to significantly reduce the uncertainty in the 
choice of the bound state potential parameters and to test the DWBA or other underlying reaction theory. In the 
combined method the ANC should be determined from an independent measurement of a peripheral reaction while 
the SF is determined from transfer reactions which are sensitive to the nuclear interior. In [20] we emphasize that 
fixing the ANC is absolutely necessary, since even when the beam energy is well above the Coulomb barrier, most of 
the reaction happens in the asymptotic region. It has been found [20] that, in (d,p) reactions, the standard single 
particle parameters for the radius r$ = 1.2 fm and diffuseness a = 0.65, which typically provide unit SF for closed 
shell nuclei, do not reproduce the value of the ANC extracted from the lower energy measurement. In this work we 
expand on the ideas of the combined method, and explore other uncertainties (such as optical potentials and higher 
order effects) to attempt a unification of the SF and the ANC, searching for a reaction description which is practical 
and gives reliable spectroscopic information. 

In section II we present a short description of the theoretical framework. In Section III we detail the results for 
our three test cases: IIIA for DWBA results using a global deuteron optical potential, IIIB for DWBA results where 
the deuteron potential is obtain from a direct fit to the elastic, IIIC for results including deuteron breakup within 
the adiabatic model, HID for checks of other higher order effects and also in HIE a further analysis of peripherality. 
Finally, a discussion of the results and conclusions are drawn in section IV. 

II. REACTION FORMALISM AND SPECTROSCOPY 

The central element of the analysis of the transfer reaction A(d,p)B [20] is the overlap function /Jt„(r) between 
bound-states of nuclei B = A + n and A which depends on r, the radius-vector connecting the center of mass of A 
with n. The square norm of the overlap function gives a model-independent definition of the SF. It is important to 
keep in mind that most of the contribution to the SF comes from the interior. 
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The radial overlap function (for B = A + n) behaves as a spherical Bessel function for large distances: 

^An(y) ( r ) r> ^ N C lj iKhi(iKr). (1) 

Here, k = y/'2 [J,An £An, where EAn is the binding energy for B — > A + n, and HAn is the reduced mass of A and 
n. R N represents a radius beyond which the nuclear potential is negligible and Cy is referred to as the asymptotic 
normalization coefficient (ANC). 

The standard practice is to take the radial dependence of the overlap function from the single particle orbital. This 
radial wavefunction is usually generated by a Woods Saxon potential with a given geometry. The depths are adjusted 
such that the single particle orbital has the correct separation energy and quantum numbers (n,l,j). Generally, the 
single particle orbital has the same asymptotic behaviour as the many-body overlap function: 

¥>An(n r ij) (r) r> « JV b nr ijiKhi(inr), (2) 

where b nT ij is the single-particle ANC (SPANC). From the relations in Eqs. ([I]) and |(2]), and under the assumption 
that the many-body overlap is indeed proportional to the single particle function all the way down to r = 0, the ANC 
and the SF are related by Cf- = S, lr ij &^ r y- So far, it is hard to check the exactness of this proportionality given the 
accuracy of the asymptotics of the wavefunction obtained from ab-initio calculations. 
For A(d,p)B, the one-step post-form finite range DWBA amplitude is given by: 

M=<^ ) /IJAy|^„^ (+) >, (3) 

where AV — V pn + V p a — U p b is the transition operator, Vy is the interaction potential between i and j, and U p b is 
the optical potential in the final-state, ip^ and -0/ ^ are the distorted waves in the initial and final states and ip pn 
is the deuteron bound-state wave function. 

In [23], the reaction amplitude is split into interior and exterior parts. The normalisation of the former is determined 

by the SF while the ANC governs the normalisation of the latter. We introduce the ratio R th (b nr ij) = a ^ pcak ^ to 

n r lj 

compare with the experimental counterpart TZ ex P = a ^i,;^ . Introducing the information about the ANC fixes the 

exterior contribution. If the reaction is completely peripheral, its cross section will scale directly as Cy, but will hold 
no information on the SF. If there is an interior contribution, the theoretical cross section has a non-trivial dependence 
on the SPANC 6„,.y which can be constrained by R th (b nr ij) = lZ exp . Hence, peripheral reactions are ideal to extract 
ANCs but useless for extracting SFs. To determine SFs one should explore non-peripheral reactions. 

The combined method presented in [13] tries to isolate the ambiguity coming from the single particle parameters 
but the resulting SFs are often smaller than those produced by shell model. As in this method, the interior part plays 
an important role to determine the SF, the results will be more sensitive to the optical potentials and coupling effects. 
With global optical potentials, the single particle parameters obtained from R th (b nr ij) — TZ exp are often far from 
the conventional values, and the corresponding angular distributions provide a worse description of the data, when 
compared to the standard procedure. The main questions we want to address here is whether DWBA (or higher order 
reaction theory) allows us to extract correct SFs when fixing the peripheral part of the transfer amplitude through 
the experimentally determined ANCs. 



III. RESULTS 

The most important consequence of the work in |'2c| is the fact that within the standard DWBA approach, the 
extracted spectroscopic factors are inconsistent with the ANCs obtained through independent measurements. This 
will be illustrated through the examples in this section. 

By standard DWBA we mean the framework in which the 1-step transfer matrix element is evaluated with incoming 
and outgoing distorted waves calculated by fitting the deuteron and proton elastic scattering with local optical 
potentials. The transfer operator contains finite range effects as well as the full complex remnant term. 

ANCs can be extracted from sub-Coulomb reactions, which are Coulomb dominated and contain virtually no 
contribution from the nuclear interior. For the 17 0(g.s.) case, we have two sets of heavy ion collision data that 
provide Cj 5/2 = 0.67 ± 0.05 fm _1 [HI] and Cj 5/2 = 0.69 ± 0.03 fm _1 [H] respectively. For the other two cases, no 

heavy ion transfer data exist, therefore alternative data was used. For 15 C(g.s.) knockout data from [23| and [13] was 
jointly used in [2|| to extract an ANC of C^ l/2 = 1.48 ± 0.18 fin -1 . Note that other values can be obtained through 
different types of reactions such as capture [2a]. We want to be consistent with the ANC extractions in all three test 
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Figure 1: 14 C(d,p) 15 C at E d = 14 MeV with a global deuteron potential (DWBAg). Data from [32 
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Figure 2: 16 0(d,p) 17 at E d = 15 MeV with a global deuteron potential (DWBAg). Data from [33 
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Figure 3: 40 Ca(d,p) 41 Ca at E d = 11 MeV with a global deuteron potential (DWBAg). Data from \34 



cases, thus we use the ANC obtained through nuclear processes. For 41 Ca(g.s.), (d,p) data at sub-Coulomb energies 
was used [131 from which we obtain C| 7 ^ 2 = 8.36 ± 0.42 fm -1 . 
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Table I: SFs and ANCs obtained from DWBA analyses: perey and perey potential for the deuteron. 





14 C (14 MeV) 


16 O (15 MeV) 


40 Ca (11 MeV) 


ro (fm) 


6 (fm -1 ) 


SF 


C 2 (fm -1 ) 


b (fm -1 ) 


SF 


C 2 (fm -1 ) 


b (fm -1 ) 


SF 


C 2 (fm" 1 ) 


1.01 


1.342 


1.40 


2.516 


0.675 


1.54 


0.700 


1.322 


1.66 


2.900 


1.11 


1.377 


1.33 


2.527 


0.753 


1.28 


0.723 


1.664 


1.18 


3.273 


1.21 


1.415 


1.27 


2.541 


0.841 


1.05 


0.745 


2.091 


0.834 


3.647 


1.31 


1.454 


1.21 


2.554 


0.940 


0.869 


0.767 


2.623 


0.583 


4.018 


1.41 


1.496 


1.15 


2.572 


1.050 


0.716 


0.789 


3.283 


0.407 


4.390 


1.51 


1.540 


1.09 


2.592 


1.173 


0.589 


0.811 


4.099 


0.284 


4.775 


1.61 


1.586 


1.04 


2.621 


1.310 


0.486 


0.833 


5.107 


0.199 


5.198 


1.71 


1.635 


0.995 


2.659 


1.462 


0.401 


0.857 


6.347 


0.141 


5.682 


1.81 


1.686 


0.953 


2.710 


1.631 


0.332 


0.883 


7.872 


0.101 


6.260 


1.91 


1.740 


0.918 


2.777 


1.818 


0.275 


0.911 


9.744 


0.073 


6.968 




r 2 


1.48 ±0.18 fm -1 


r 2 

^ exp 


0.67 ±0.05 fm -1 


r 2 

v - y exp 


8.36 ± 0.42 fm -1 



The transfer calculations are performed for 14 C(d,p) 15 C(g.s.) at E l d ab = 14 MeV [33], 16 0(d,p) 17 0(g.s.) at E l d ab = 15 
MeV [331 , and 40 Ca(d,p) 41 Ca(g.s.) at E l d ab — 11 MeV [3J|. These reactions were chosen because both elastic and 
transfer data exist at the same energy. All data can be found on the database website at Michigan State University 
[35|]. Calculations were performed with the code FRESCO (36| . 



A. DWBA with global deuteron potential 



We perform finite-range post-form DWBA calculations, including the full complex transition operator AV. The 
Reid-Soft-Core interaction [30|] is used for the deuteron ground state wavefunction, and includes both S- and D- waves. 
We use the central part of the Reid-Soft-Core interaction for V np in the transfer operator. The diffuseness of the single 
particle orbital of the final state is kept fixed (a=0.65 fm), but the radius is varied to generate a range of SPANCs. 
Throughout this work we will always use the CH89 global parameterization for the outgoing distorted waves. 
We have checked that by using a different proton global potential, spectroscopic factors change less than 10 %. We 
will look into several approaches to determine the initial wavefunction. In this subsection we take the Perey & Perey 
deuteron potential (PP) [HI]. 

The resulting angular distributions are displayed in Figs. [U [2] and 02 for a subset of ro, and compared to data. 
We will refer to these calculation as DWBAg. Even though it provides a fair description of the first peak of the 
distribution, it is inadequate for the large angles, where higher order effects become important. This is a well known 
characteristic of DWBA. 

For each single particle radius ro of B = A ± n, the normalization of the first transfer peak is used to determine 
the spectroscopic factor, and the ANCs are obtained directly from Cf- — S Ur ij b^ij- Results are presented in Table 
HI If one takes the standard radius of 1.21 fm, one can see that the extracted SFs are all close to unity, but none of 
the ANCs are consistent with the values extracted from experiment, with the most serious mismatch for 41 Ca, nearly 
a factor of 2. Knockout measurements suggest that in 17 and 41 Ca there is a reduction of the SF to ss 0.6. One 
can see from Table IU that this SF is reproduced at r w 1.4 fm for 17 and r ~ 1.3 fm for 41 Ca, but the SF/ANC 
inconsistency is not resolved. 

Finite range effects are known to be important in (d,p) reactions, but remnant contributions are often assumed 
to be small. We have repeated DWBAg calculations excluding the remnant term and find that indeed, remnant 
contributions are insignificant for the reactions of 14 C and 16 0, but can contribute up to 30% for the 40 Ca(d,p) 
reaction. 



B. DWBA with fitted deuteron potential 



The elastic channel data is also available for the three reactions under scrutiny, 14 C(d,p) 15 C(g.s.) at E l d ah = 14 
MeV, 16 0(d,p) 17 0(g.s.) at E l d ab = 15 MeV and 40 Ca(d,p) 41 Ca(g.s.) at E l d ab = 11 MeV. We have repeated this study 
using optical potentials fitted specifically to the corresponding elastic data, using SFRESCO [36]. The optical potentials 
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Table II: Deuteron optical potential parameters resulting from the fit to elastic data (fit 1). 



A 


E beam (MeV) 


V R {MeV) 


Rr (fm) 


an (fm) 


W d {MeV) 


Ri (fm) 


ai (fm) 


K (MeV) 


Rso (fm) 


a so (fm) 


14 


14 


92.880 


1.1486 


0.596 


3.848 


1.085 


1.347 


9.164 


0.883 


0.244 


16 


15 


115.911 


1.017 


0.846 


11.257 


1.073 


0.584 


11.600 


0.578 


0.343 


40 


11 


115.177 


1.040 


0.712 


5.287 


1.375 


0.856 


4.486 


0.382 


0.266 



Table III: SFs and ANCs obtain from DWBA analyses: deuteron optical potential directly fitted to the corresponding elastic 
data. 





14 C (14 MeV) 


16 (15 MeV) 


40 Ca (11 MeV) 


ro (fm) 


b (fn^ 1 ) 


SF 


C 2 (fm" 1 ) 


b (fm" 1 ) 


SF 


C 2 (fm" 1 ) 


b (fm" 1 ) 


SF 


C 2 (fm- 1 ) 


1.01 


1.34 


1.50 


2.70 


0.67 


1.38 


0.632 


1.32 


1.23 


2.159 


1.11 


1.37 


1.42 


2.70 


0.75 


1.09 


0.622 


1.66 


0.985 


2.729 


1.21 


1.41 


1.35 


2.71 


0.84 


0.869 


0.616 


2.09 


0.748 


3.272 


1.31 


1.45 


1.27 


2.70 


0.93 


0.693 


0.613 


2.62 


0.545 


3.754 


1.41 


1.49 


1.20 


2.69 


1.04 


0.555 


0.613 


3.28 


0.387 


4.172 


1.51 


1.54 


1.13 


2.68 


1.17 


0.447 


0.617 


4.10 


0.271 


4.557 


1.61 


1.58 


1.05 


2.66 


1.31 


0.363 


0.625 


5.11 


0.189 


4.945 


1.71 


1.63 


0.988 


2.64 


1.46 


0.297 


0.636 


6.35 


0.133 


5.380 


1.81 


1.68 


0.921 


2.62 


1.63 


0.244 


0.651 


7.87 


0.095 


5.904 


1.91 


1.74 


0.859 


2.60 


1.81 


0.202 


0.671 


9.75 


0.069 


6.555 




r 2 


1.48 ±0.18 fm -1 


r 2 

^ exp 


0.67 ± 0.05 fm" 1 


r 2 

l ~ y exp 


8.36 ± 0.42 fm -1 



resulting from the elastic fits are presented in Table ITT1 and differ somewhat from the global compilations. The resulting 
transfer cross sections, together with the corresponding elastic scattering, are shown in Figures |4j [5] and [6l We will 
refer to these calculations as DWBAf. For both the 14 C(d,p) and the 41 Ca(d,p) reactions, the SF/ANC inconsistency 
remains whereas for the 16 0(d,p) a spectroscopic factor of 0.6-1 produces ANCs much closer to the value extracted 
from a sub-Coulomb measurement. This can be seen from Table UlTl where SPANCs, SFs and ANCs are presented as 
a function of the single particle radius for the three reactions under study. 

C. Adiabatic deuteron potential 

It is well known that the deuteron breakup can have an effect on the transfer cross sections. One way to take 
this into account would be to couple explicitly the deuteron continuum within the Continuum Discretized Coupled 
Channel (CDCC) method [3§|. The complexity of the procedure would not be practical for the general experimental 
community. Alternatively, the adiabatic method developed by Johnson and Soper [2§| can be used to describe the 
entrance channel and obtain the transfer cross section while including the deuteron breakup channel. This procedure 
is by far simpler than the CDCC approach and can easily be used in systematic studies (e.g. [lfl The three 

body wavefunction for the d — T system obtained with the adiabatic approximation has the correct properties in the 
range of V np yet should not be used in the asymptotic region. For this reason we perform finite range calculations but 
do not include the remnant term in Eq.j3]). Also important is that the deuteron wavefunction be an eigenstate of V np 
used in the transfer operator. Transfer matrix elements with the full RSC potential are not simple. Thus we have 
used the central gaussian, refitted to reproduce the correct binding energy and the same D as the RSC potential, for 
calculating both, the wavefunction and the transfer operator. 

We use the CH89 nucleon potential for the U n r and U p t at half the deuteron incident energy and calculate the 
adiabatic potential using the parametrization in [4(j] which included finite range corrections. We will refer to these 
results as the adiabatic wave approximation (ADWA) [4lJ. The cross section obtained for the three reactions under 
study are presented in Figs. [7][9l In table IIVI we show the results for the extracted SF and corresponding ANCs. 
An overall reduction of the spectroscopic factors is observed. The angular distributions show some improvement as 
compared to those within DWBAg or DWBAf. However, it becomes clear that deuteron breakup is unable to remove 
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Table IV: SFs and ANCs obtained from ADWA analyses. 





14 C (14 MeV) 


16 (15 MeV) 


40 Ca (11 MeV) 


ro (fm) 


6 (fm -1 ) 


SF 


C 2 (fm -1 ) 


b (fm -1 ) 


SF 


C 2 (fm -1 ) 


6 (fm -1 ) 


SF 


C 2 (fm" 1 ) 


1.01 


1.34 


1.16 


2.08 


0.67 


1.45 


0.659 


1.32 


1.80 


3.14 


1.11 


1.37 


1.11 


2.10 


0.75 


1.19 


0.675 


1.66 


1.30 


3.57 


1.21 


1.41 


1.06 


2.12 


0.84 


0.974 


0.688 


2.09 


0.895 


3.92 


1.31 


1.45 


1.01 


2.14 


0.93 


0.794 


0.701 


2.62 


0.609 


4.19 


1.41 


1.49 


0.969 


2.17 


1.04 


0.646 


0.712 


3.28 


0.410 


4.43 


1.51 


1.54 


0.928 


2.20 


1.17 


0.526 


0.724 


4.10 


0.278 


4.67 


1.61 


1.58 


0.891 


2.24 


1.31 


0.429 


0.737 


5.11 


0.189 


4.94 


1.71 


1.63 


0.858 


2.30 


1.46 


0.351 


0.751 


6.35 


0.131 


5.28 


1.81 


1.68 


0.830 


2.36 


1.63 


0.289 


0.768 


7.87 


0.092 


5.72 


1.91 


1.74 


0.807 


2.44 


1.81 


0.238 


0.788 


9.75 


0.066 


6.30 




r 2 


1.48 ±0.18 fm -1 


r 2 


0.67 ±0.05 fm -1 


r 2 

v - y exp 


8.36 ± 0.42 fm -1 



the inconsistency between SF and ANC, specially in the Ca case, where a persistent factor of 2 remains. 

As mentioned above, the ADWA results here presented do not include the remnant. Results from DWBAg show 
that remnant contributions are not important for the lighter cases under study. Even though remnant contributions 
to 40 Ca(d,p) 41 Ca are not negligible, the magnitude is much smaller than the mismatch observed and we do not 
expect their inclusion within ADWA would change the conclusions. Also, the deuteron wavefunction for the adiabatic 
calculations is s-wave only (does not contain the d-wave and the tensor interaction). This effect alone reduces the 
spectroscopic factors by ~ 10%. However our conclusions, namely that deuteron breakup cannot account for the 
inconsistency between SF and ANC, remains. 

D. Other higher order effects 

The transfer couplings for the examples we are studying are relatively strong, therefore one solution to the incon- 
sistency between SF and ANC could reside in higher order processes in the reaction mechanism, other than deuteron 
breakup. In these calculations our starting point is DWBAf. First we consider multiple transfer couplings within the 
Coupled Reaction Channel approach (CRC). Note that some of these effects are accounted for within ADWA. We 
restrict ourselves to the (d,p) transfer coupling connecting ground states of the A and B. For all three reactions, we 
increase the number of iterations until convergence is achieved. We refit the deuteron optical potentials, so that the 
elastic scattering is still well reproduced. 

• For 14 C(d,p) 15 C(g.s.) at E l ^ b = 14 MeV, we find that there is a significant effect of higher order transfer 
couplings in the cross section, which amounts to a reduction of the required SF. Still this does not completely 
solve the inconsistency. For example, for ro = 1.3 fm, S — 1.013 and C 2 = 2.14 fm -1 . 

• CRC effects for 16 0(d,p) 17 at 15 MeV are weak (a few percent). 

a p or 40 Ca(d,p) 41 Ca at 11 MeV, CRC increases the transfer cross section by ~ 20% which produces lower SFs. 
This makes the inconsistency SF/ANC more severe. 

Another source of higher order effects comes from explicit multi-step excitations of the target. These targets are 
spherical but vibrate, and thus there are several possible couplings that could be considered within the coupled channel 
Born approximation (CCBA). Taking into account the nuclear deformation lengths determined from nucleon inelastic 
studies |42l. |43|. we have included the excited state of the target that couples most strongly to the ground state, the 
3 _ state. Again the deuteron optical potential is refitted to reproduce the elastic scattering correctly. We find that 
for 14 C, these couplings have very little effect. This is not the case for both 16 O and 40 Ca. 

Let us first consider 17 0. Its ground state is described in terms of a coupled-channel equation that generates a 
dominant c? 5 / 2 component. Changing ro can have a dramatic effect in the structure composition (mixing in p- and 
f- waves), and thus one needs to limit the range of single particle parameters where a realistic structure is preserved. 
Results are shown in table |V] Including the coupling to the 3 _ state, reduces the transfer cross section, increasing 
the SF. The resulting ANC of the d 5/2 for a unit SF, is 0.56 fm" 1 . 
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Table V: Variation of the SFs and the ANCs with r for CCBA calculations of 16 0(d,p) 17 at E d = 15 MeV. 



rO 


SF 


b(ld5/2) 


C 2 (d5/2) 


1.16 


1.28 


0.6795 


0.59 


1.20 


1.17 


0.7003 


0.58 


1.24 


1.08 


0.7215 


0.56 


1.28 


0.99 


0.7430 


0.55 


1.32 


0.92 


0.7647 


0.53 



Table VI: Variation of the SFs and the ANCs with r for CCBA calculations of 40 Ca(d,p) 41 Ca at E d = 11 MeV. 



rO 


SF 


b(lf7/2) 


C7 2 (/7/2) 


1.22 


1.16 


1.7869 


2.07 


1.24 


1.09 


1.8541 


3.75 


1.26 


1.03 


1.9233 


3.81 


1.28 


0.97 


1.9944 


3.85 


1.30 


0.91 


2.0674 


3.88 



In the Ca case, the situation is different. Taking again only values of ro that generate coupled-channel wavefunc- 
tions for 41 Ca that are still dominantly /V/2, one obtains transfer cross sections that are smaller than in the single 
particle case. Thus the extracted SFs are higher (see table IVI|) . Nevertheless, the ANCs associated with the /V/ 2 
component are still much smaller than the number extracted from the sub-Coulomb measurement. We have performed 
CRC iterations on CCBA calculations for 40 Ca(d, p) 41 Ca at E d = 11 MeV. In this case the transfer couplings involved 
are: d+ 4Q Ca(g.s.) — > p+ 41 C&(g.s.) and d+ 40 Ca(3^) — > p+ 41 Ca,(g.s.), with spectroscopic factors as in table IVll We 
find the picture does not change significantly. 

E. DBWA with a simultaneous fit 



Table VII: Single particle orbitals adjusted to produce the experimental ANC, when SF is unity. 





15 C 


17 


41 Ca 


ro (fm) 


1.10 


1.25 


1.30 


a (fm) 


0.52 


0.60 


0.75 


C 2 (fm" 1 ) 


1.62 


0.68 


8.24 


Rrmsifm) 


5.1 


3.5 


4.2 



It is impossible to rule out other excitation mechanisms which may contribute to the transfer process, however we 
do not expect these will be stronger than the ones discussed in section HID. One interesting question is whether there 
would be any way of obtaining the desired consistency between SF and ANC within the simplified DWBA picture, 
specially in the 41 Ca case, where the mismatch is so large. We assume that we can introduce all higher order processes 
into a local effective deuteron optical potential. We impose that the SF for a neutron outside a closed shell be unity, 
and chose the single particle parameters such that they reproduce the ANC extracted from experiment. We perform 
a nine parameter fit of the deuteron optical potential to both the elastic and the transfer data. Our starting potential 
is that obtained in section UlIBl The fit consists of a standard \ 2 minimization procedure and the code SFRESCO is 
used [36]. For each case we do find a potential that is able to reproduce the elastic and the transfer simultaneously, 
under the constraint of SFs and ANCs. The resulting real (Vfit2) and imaginary (Wfit2) potentials from this fit are 
presented for all three reactions, as the dashed lines (fit2) in Figs. [10l The parameters are explicitly shown in table 
I Villi The solid lines correspond to the previous elastic fit (fitl), whose parameters were shown in table [ill It appears 
that for 14 C + d at E d = 14 MeV, the real part becomes slightly more diffuse, whereas for 41 Ca + d at E d = 11 MeV, it 
becomes clearly less diffuse. The 16 + d at E d = 15 MeV case is rather unchanged, which reflects the fact that to start 
with the inconsistency in this case was minor. We find that the most significant change in the resulting potentials is 
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Table VIII: Deuteron optical potential parameters resulting from the simultaneous fit to elastic and transfer data. 



A 


E beam (MeV) 


Vk(MeV) 


Rr (fm) 


a,R (fm) 


Wd(MeV) 


Ri (fm) 


ai (fm) 


l4o(MeV) 


Rso (fm) 


aso (fm) 


14 


14 


113.015 


0.902 


0.875 


10.659 


1.007 


0.900 


14.342 


1.281 


0.694 


16 


15 


111.531 


0.948 


0.871 


10.211 


1.635 


0.437 


10.991 


0.505 


0.541 


40 


11 


126.922 


1.004 


0.500 


5.625 


1.800 


0.736 


6.606 


0.988 


0.4038 



an increase of the imaginary part in the surface region, which is obtained either by a large imaginary depth or by a 
shift toward the surface. This suggests that the standard DWBA is overestimating the surface contribution. 

Changes in the deep interior do not reflect sensitivity of the transfer process to this region. Transfer reactions are 
known to be surface peaked. This usually means that impact parameters smaller than the sum of the projectile and 
target radii do not contribute to the reaction. In A(d,p)B, this is related to the cutoff on RdA- The results of this 
standard peripherality check is shown in Fig. [TT]for the three reactions under study. The percentage ratio of the cross 
section integrated up to RdA — R to the total transfer cross section R(%) is plotted. In the deep interior this ratio is 
zero and it goes to 100 % for large distances. The region where it increases rapidly corresponds to the surface region, 
where the transfer takes place. 

However, to probe the sensitivity to the SF, one needs to analyse the dependence with R n A- Under a zero range 
approximation of the deuteron, these two tests would be identical. When taking into account the finite range of V np 
and the remnant part of the full transition operator, peripherality in RdA is not equivalent to peripherality in R n A- 
To illustrate this fact, peripherality tests were performed by evaluating the interior contribution to the total transfer 
cross section relative to R n A- These peripherality tests are based on DWBAf (fit 1) and DWBA with the deuteron 
potential fit simulaneously to elastic and transfer (fit 2). We take the radial integrals in the coordinates of R n A, RpB, 
and truncate the integration in R n A to a maximum value r: 

i poo pr 

^ (r) 7o Jo ^InA^^dRnAdRBp) 2 . (4) 

The percentage ratio of this value to the full integration is presented in Fig. [TJl As before, if r = fm, the ratio 
should be zero, and it should tail off at 100 % when r becomes very large. Results show that, in all cases there is 
no contribution up to 2 fm. One can also see that, even though 14 MeV is above the Coulomb barrier, the 14 C(d,p) 
happens at rather large distances, due to its loosely bound nature. Contributions up to 20 % from the surface/interior 
are present in both 16 + d at E d = 15 MeV and 41 Ca + d at E d = 11 MeV. 

IV. DISCUSSION AND CONCLUSIONS 

In this work we have studied transfer reactions to states considered good single particles, with three different De- 
values, namely 14 C(d,p) 15 C(g.s.) at E l d ab = 14 MeV, 16 0(d,p) 17 0(g.s.) at E l f> = 15 MeV 40 Ca(d,p) 41 Ca(g.s.) at 
E l d ab = 11 MeV. All these reactions are above the Coulomb barrier and therefore contain some information from the 
interior. The standard DWBA method, using global optical potentials and the typical single particle parameters, 
produces SFs in agreement with shell model predictions, however the corresponding ANCs are not consistent with 
those extracted from independent measurements. If one imposes, within the DWBA formulation, ANCs that are 
consistent with the experimental values, the extracted SFs are significantly reduced compared to the shell model 
predictions. 

Some improvements on the SF/ANC mismatch can be obtained by using a deuteron optical potential fitted directly 
to the corresponding elastic data, at the relevant energy. In particular, for 16 0(d,p) we obtain SF/ANC consistency. 
However the problem for the other two cases is not resolved. The deuteron adiabatic potential, which takes into 
account breakup, can change the SF up to 30 %. This improves the situation for 14 C(d,p) but fails to bring the 41 Ca 
ANC anywhere close to the value extracted from an independent sub-Coulomb measurement. 

Deuteron breakup is not the only important ingredient to the solution of this problem: in all cases we find an 
influence of surface higher order effects. For 14 C(d,p) 15 C(g.s.) at E l d ab — 14 MeV, CCBA effects are negligible but 
CRC effects produce a reduction of the ANC consistent with unity SF: C 2 — 2.14 fm" 1 , a significant reduction 
compared to the DWBA case, which together with the reduction from deuteron breakup produces an ANC much 
closer to the experimental value C 2 = 1.48 ±0.18 fin" 1 . For 16 0(d,p) 17 0(g.s.) at E l d ab = 15 MeV, CRC effects are 
small but CCBA couplings to the 3~ state in 16 were found to be important. This brings down the C 2 from 0.76 to 
0.56, neither too far from the experimental number 0.67 ± 0.05 fm -1 . Finally, the 40 Ca(d,p) 41 Ca(g.s.) case remains 
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problematic: both CRC and CCBA are relevant but act in the wrong direction. Overall, an ANC, consistent with 
unity SF, falls short by nearly a factor of 2. 

Contrary to (e,e'p) measurements, transfer reactions are surface peaked and it is disconcerting that the traditional 
methods to handle higher order effects at the surface are not able to solve the SF / ANC discrepancy for one of our 
test cases. The very fact that, even when the energies are well above the Coulomb barrier, there is such a large 
contribution from the peripheral region, makes it extremely important to pin down the ANC input unambiguously. 
We cannot rule out the possibility of a problem in the 41 Ca ANC we extracted from other data. Experiments to 
measure ANCs for this case is crucial to settle this matter. In the future we suggest that experiments be designed for 
the extraction of ANCs in parallel with the corresponding experiments aimed at extracting SFs. 
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Figure 4: a) Elastic scattering fit 14 C + d; b) 14 C(d,p) 15 C at E d = 14 MeV with a fitted deuteron potentials (DWBAf). Both 
elastic and transfer data from [321 ]. 





e cm (deg) e cm (deg) 



Figure 6: a) Elastic scattering fit 14 C + d; b) 40 Ca(d,p) 41 Ca at E d = 11 MeV with a fitted deuteron potentials (DWBAf). 
Elastic data from [H] and transfer data from (34)]. 
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Figure 7: 14 C(d,p) 15 C at E d = 14 MeV with ADWA. Data from [12]. 




Figure 8: 16 0(d,p) 17 at E d = 15 MeV with ADWA. Data from 




Figure 9: 40 Ca(d,_p) 41 Ca at E d = 11 MeV with ADWA. Data from (H. 
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Figure 10: Optical potentials: fits to elastic scattering only - real (Vfitl) and imaginary (Wfitl) - and simultaneous fits to the 
elastic and transfer channel - real (Vfit2) and imaginary (Wfit2). a) 14 C + d at E d = 14 MeV, b) 16 + d at E d = 15 MeV, c) 
41 Ca + d at E d = 11 MeV. 




Figure 11: Peripherality test for the reactions (using a cut in the radial distance between projectile and target): a) 14 C + d at 
E d = 14 MeV, b) 16 + d at E d = 15 MeV, c) 41 Ca + d at E d = 11 MeV. 
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Figure 12: Peripherally test for the reactions (using cutoff in the radial distance R n A). a) C + d at Ed = 14 MeV, b) O + d 
at E d = 15 MeV, c) 41 Ca + d at E d = 11 MeV. 



